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ABSTRACT 

We present 0.55 x 10 6 -particle simulations of the accretion of high-density dwarf galaxies by low- 
density giant galaxies, using models that contain both power-law central density cusps and point masses 
representing supermassive black holes. The cusp of the dwarf galaxy is disrupted during the merger, 
producing a remnant with a central density that is only slightly higher than that of the giant galaxy 
initially. Removing the black hole from the giant galaxy allows the dwarf galaxy to remain intact and 
leads to a remnant with a high central density, contrary to what is observed. Our results support the 
hypothesis that the persistence of low-density cores in giant galaxies is a consequence of supermassive 
black holes. 

Subject headings: galaxies: evolution — galaxies: elliptical and lenticular, cD galaxies: interactions — 
galaxies: nuclei 



1. INTRODUCTION 

Bright spheroids (elliptical galaxies and the bulges of 
spiral galaxies) tend to be less dense than faint spheroids, 
in two distinct ways. Bright spheroids have larger scale 
lengths, r e ~ L 3 / 4 , and lower mean densities within r e , 
where r e is the half-light radius and L is the total lumi- 
nosity. Bright galaxies also tend to have density profiles 
that are less concentrated toward the center. The luminos- 
ity densities of elliptical galaxies and bulges rise approx- 
imately as power laws at the smallest observable radii, 
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Faint galaxies (My > -20) have 1.5 < 7 < 2.5 



while bright galaxies have < 7 < 1.5 (|Crane et al. 1993 



Fcrrarese et al. 1994 |Gebhardt ct al. 1996[ ) 



The existence of the second correlation raises two ques- 
tions. Why do bright galaxies have lower central concen- 
trations than faint galaxies; and how is this correlation 
maintained in the face of mergers? Big galaxies will some- 
times accrete small galaxies, and the small galaxy wil l be 
resistant to ti d al disruption due to its h igher density ( |Ko-| 
rmendy 1984 ; Balcells fc Q uinn 1989[ ). A small galaxy 
that survived such a merger with its central regions intact 
would create a new, high-density core in the merger rem- 
nant, destroying the observed correlation between L and 
7- 

The survivability of compact galaxies during mergers is 
relevant to the origin of galaxies with kinematically dis- 
tinct cores. Forbes, Franx & Illingworth (1995) investi- 
gated whether the properties of such galaxies were consis- 
tent with the accretion hypothesis. They found no strong- 
evidence for abnormal surface brightness profiles ('cores 
within cores') but noted that "This problem can possi- 
bly be circumvented by adding a massive black hole to the 
host galaxy. The [host] galaxy might be able to disrupt the 
victim galaxy, resulting in a density profile more like that 
which is seen." They suggested that "Simulations with 
realistic density profiles and possible black holes would be 
valuable." 

Here we present the results of such simulations. We 
consider mergers between initially spherical galaxies with 



a mass ratio of 10 : 1; each galaxy has a power-law central 
density cusp as well as a central point mass representing 
a supermassive black hole. The galaxies' structural pa- 
rameters are scaled in accordance with observed relations 
(§2). We find, in agreement with the suggestion of Forbes, 
Franx & Illingworth (1995), that the black hole in the 
primary galaxy is effective at disrupting the cusp of the 
secondary galaxy, yielding a remnant density profile that 
is only slightly more centrally concentrated than that of 
the primary galaxy at the start of the simulation. 

2. METHOD 

The initial galaxies were generated from 
Dehnen's (1993) spherical density law, p(r) = 
[(3 - 7)M/47ra 3 ] (r/a)"T(l + r/a)~<~ 4 . Initial particle 
velocities were assigned from the unique isotropic distri- 
bution function that reproduces Dehnen's density law in 
the combined potential of the st ars and a central poin t 
mass representing the black hole (Iremaine ct al. 1994). 
Our primary galaxy had 7 = 1, characteristic of the shal- 
low cusps of bright galaxies, and our secondary galaxy had 
7 = 2, characteristic of the steep cusps of dwarf galaxies. 
Each "black hole" was assigned a mass of 2 x 10~ 3 times 
that of its parent galaxy. This is slightly above, but con- 
sistent with, the current best estimate of ~ 0.0012 for the 
average ratio of bla ck hole mass to luminous g alaxy mass in 
the local universe ( |Merritt fc Fcrrarese 2001 ). Henceforth 
the subscripts 1 and 2 will refer to the primary (massive) 
and secondary (dwarf) galaxies respectively. All of the 
simulations had M1/M2 — 10. 

The mass M\ and scale length a\ of the primary galaxy 
were set to unity. The orbital period at the half-mass ra- 
dius of the primary is ~ 33.3 in model units. The scale 
length of the secondary galaxy was computed from the 
relation 
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where the ratio r e /a between effective radius and scale 
length is given by Dehnen (1993) as a function of 7. Real 
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spheroids have r e oc L a and M/L oc L* with M/L the 
mass-to-light ratio of the stars; hence r P oc M a /(^ +1 ), Set- 



ting a » 0.75 (e.g. |Valluri fc Merritt 1998|) and /3 w 0.25 
(e.g. |Faber et al 1987ft gives r e , 2 /r e ,i ?a (M2/M1) 3 / 5 . 
Using our assumed mass ratio of 10 : 1 and the adopted 
values for 7, equation (1) then implies 02/01 = 0.618. 

The center of the secondary galaxy was initially dis- 
placed from the center of the primary galaxy by ~ 3 
times the half-mass radius of the primary along the y- 
axis. The initial velocity of the secondary galaxy was 
assigned one of four values determined by k = L/L C i rc , 
the angular momentum of the secondary's orbit in units 
of the angular momentum of a circular orbit. We chose 
k = {0,0.2,0.5,0.8}; the initial velocity vector of the sec- 
ondary galaxy was parallel to the x axis. 

The primary galaxy was assigned N\ — 0.5 x 10 6 parti- 
cles and the secondary galaxy N2 — 0.05 x 10 6 particles. 
Since N1/N2 = M\/M% — 10, all particles had the same 
mass. This was done to avoid any spurious relaxation due 
to mass segregation. 

The evolution was followed using the tree code GADGET 
of Springcl ct. al (2000), a parallel algorithm with continu- 
osly variable time steps for each particle. The parameters 
in this code that most strongly affect the accuracy of the 
integrations are the softening length h and the time step 
parameter r\. We chose h = 8 x 10~ 4 , slightly smaller than 
the radius ~ 10~ 3 at which the two black holes would 
be expected to form a hard binary. This softening length 
corresponds to < 1 pc with typical scalings. Our value 
?/ = 0.008 for the accuracy parameter was chosen so as to 
require > 10 2 time steps for a star in a circular orbit with 
radius h around the large black hole. Extensive testing 
was carried out to ensure that the selected values for the 
tree-code parameters resulted in no changes in the density 
profiles of either galaxy on length scales > 10~ 3 when inte- 
grated in isolation. The merger simulations were continued 
until roughly one crossing time of the primary galaxy after 
formation of the hard black-hole binary. Integrations were 
carried out on the parallel supercomputers at the Rutgers 
Center for Advanced Information Processing (a Sun HPC- 
10000) and at the San Diego Supercomputer Center (Cray 
T3E). All simulations used 16 processors. 

In order to test the influence of the secondary's black 
hole on the evolution, each run was repeated using a sec- 
ondary galaxy with the same initial density profile but 
lacking a central point mass. To test the importance of 
the primary's black hole, a further set of runs were car- 
ried out in which neither galaxy contained a black hole. 
The velocity distribution functions for the models without 
black holes were computed as in Dehnen (1993). 

3. RESULTS 

The results are summarized in Figures 1 and 2. Figure 
la shows, for k — 0.2, the number of particles initially as- 
sociated with the two galaxies that remain within spheres 
of ra dius 0.1 (primary') and 0.05 (secondary) around their 



respe ctive black holes; these are roughly the respective 
radii inside of which the final density profile of each galaxy 
shows substantial evolution. Figure 1 reveals two ways in 
which the merger causes the central density of the final 
object to be less than the sum of the central densities of 
the two initial objects. 1. The density of stars around the 



black hole in the primary galaxy drops steadily through- 
out the simulation, and drops suddenly whenever the sec- 
ondary galaxy passes near the primary's center; the drops 
are preceded by jumps resulting from the impulsive in- 
crease in the gravitational force caused by the passage of 
the small galaxy. 2. The density of stars around the black 
hole in the secondary galaxy also drops, but only after the 
secondary passes sufficiently close to the center of the pri- 
mary. In the case of plunging orbits, k = or 0.2, the 
first close passage occurs during the first infall; the drop 
in density of the secondary is again preceded by a jump 
indicating that the mechanism is impulsive addition of en- 
ergy, this time from the gravitational force of the primary's 
black hole. In the case of the more circular orbits, n = 0.5 
and 0.8, the central density of the secondary remains high 
until the radius of its orbit has decayed to ~ 0.1. 

The presence of the black hole in the secondary galaxy 
might naively be expected to contribute to the stability 
of that galaxy's core by enhancing the gravitational force 
there. In fact the opposite is true: the density at late 
times in the runs without a second black hole is generally 
higher than in the run with two black holes (Fig. la). 
This is probably a consequence of three-body ejection of 



stars by the pair of black holes (Quinlan 1996). However 
the softening of the inter-partic e forces in GADGET does 
not permit this process to be followed with high accuracy. 
The final central densities that we observe in the runs with 
two black holes should therefore be interpreted as upper 
limits. 

Removing the black hole from the primary galaxy has 
almost no effect on its evolution, but the central density 
of the secondary galaxy now remains high throughout the 
merger (Figure lb). 

Figures 2a-d show the final density profiles from the four 
runs in which both galaxies contained black holes. The 
profiles are centered on the black ho le binary and were 
computed using MAPEL ( Merritt 1994 ). The secondary's 
density (Fig. 2a) is lowered at small radii, r < 0.05, due 
to impulsive heating from the primary black hole, and at 
large radii, r > 0.1, due to tidal forces from the primary 
galaxy. The greatest decrease in central density occurs 
in the runs with lowest k, i.e. the most plunging orbits. 
In the case of stars initially associated with the primary 
galaxy, the final density profile is remarkably independent 
of k (Fig. 2b). The primary galaxy was found to always 
exhibit a p ~ r -o.5 cen tral density cusp at the end of 
the simulations; a similar result was found by Nakano & 
Makino (1999a, b) when they dropped a black hole into the 
core of a galaxy containing no black hole. The net result 
(Fig. 2c) is a remnant density profile that rises modestly 
above that of the initial, primary galaxy at small radii, 
more so in the runs with larger k. However the slope of 
the surface brightness profile (Fig. 2d) is only slightly 
increased near the center compared with that of the pri- 
mary galaxy. We note that many bright elliptical galaxies 
also show mild inflections in their central surf ace bright- 
ness profiles o n similar scales (e.g. Figure 1 of Merritt & 



Fridman 1995), perhaps relics of past mergers. 



When the black hole is removed from the primary 
galaxy, the secondary galaxy survives the merger intact, 
with almost no change in its density profile at small radii 
(Fig. 2e). The final density profile of the remnant is 
roughly a superposition of the two initial profiles (Fig- 
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ure 2g); a substantial inflection in the surface brightness 
profile now appears at a radius of ~ 0.1 (Figure 2h). Pro- 
files of this sort ar e rarely if ever observe d in the brightest 
elliptical galaxies (Ferrarese et al. 1994; Gcbhardt et al 
19960 



4. DISCUSSION 



Our merger simulations are the first in which the cen- 
tral d ensities of the merging galaxies are realistically high 
and ii | which both galaxies contain "black holes." Barnes 



conclusion that all but the most plunging orbits resulted 
in merger remnants with unphysically steep cusps. Holley- 
Bockclmann & Richstone also observed a stronger depen- 
dence of the secondary's density evolution on its initial 
orbit than seen here, probably because of the approximate 
way in which orbital decay was treated. 

Our results support the hypothesis that the persistence 
of low-density cores in giant galaxies is a consequence of 
the existence of sup ermassive black holes ( Forbes, Franx 
& Illingworth 1995). The central density did neverthc- 



(1999) described mergers of identical Dchncn models with 
cusp slopes 7 = {1,1.5,2} but without black holes; he 
found that the cusps survived the mergers with little 
change. We find a similar result here in the runs without 
black holes except that the density of the primary galaxy 
falls due to energy input from the secondary. The oppo- 
site case - mergers between galaxies with central black 
holes but no cusps - has been treated extensively (e.g. 



Ebisuzaki, Makino fc Okumura 1991 ; Governato, Colpi & 
Maraschi 1994| ; [Making 1997| ). A number of authors have 



investigated the effects of dropping one or two black holes 
into a pre-existing galaxy, either with a cent ral density 
cusp (Quinlan fc Hcrnquist 1997 ) or without ( Nakano ~&| 



Makino 1999a, 1999b). The simulations closest in spirit to 
ours are those of Holley-Bockelmann & Richstone (2000). 
These authors followed the evolution of a dense secondary 
galaxy as it moved in the fixed gravitational field of a less- 
dense primary containing a black hole; the secondary con- 
tained no black hole. The decay of the secondary's orbit 
was induced artificially using standard expressions for the 
dynamical friction, with corrections due to the changing 
mass of the secondary. Since the density of the primary 
galaxy was fixed in these simulations, the authors did not 
observe the drop in the primary's central density which 
we see here (Fig. 2b) ; hence the central densities of their 
merger remnants were too high, partly explaining their 



less increase in our simulations, and since a typical bright 
galaxy is expected to hav e accreted many sma ller galaxies 
since its formation (e.g. Lacey & Cole 1993), one might 
still predict the formation of dense nuclei in bright galax- 
ies, contrary to what is observed. However our simulations 
are not able to follow the evolution of a black-hole binary 
on sub-parsec scales as it ejects star s from the co re and 
lowers the density there still more ( |Quinlan 199^ ). Fu- 
ture work should be directed toward understanding the 
effects of hierarchical mergers on galaxy density profiles, 
using TV-body codes that can deal efficiently with unsoft- 
ened particles and hence follow the interaction of a binary 
black hole with the surrounding stars. 
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Fig. 1. — Time evolution of the merger simulations with k = 0.2. (a) Runs with a black hole in the primary galaxy. Top panel: separation 
between the primary and secondary black holes, in x (thick line), y (intermediate line), and z (thin line). Middle panel: number of particles 
initially associated with the primary galaxy in a sphere of radius 0.1 centered on its black hole. Bottom panel: filled circles indicate number 
of particles initially associated with the secondary galaxy in a sphere of radius 0.05 centered on its black hole. Open circles correspond to 
runs in which the secondary galaxy contained no black hole; N2 is the number of particles in a sphere of radius 0.05 centered on the point 
of maximum density in the secondary galaxy, (b) Similar to (a) except that neither galaxy contained a black hole. Centers of both galaxies 
were defined as points of maximum density. 
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Fig. 2. — Final density profiles from the simulations in which each galaxy contained a black hole (a-d) and in which neither galaxy contained 
a black hole (e-g). The four thin curves in each frame correspond to the four different orbits; the thin curve extending farthest to the left 
corresponds to re = 0.8. (a), (e) Space density of stars initially associated with the secondary galaxy; thick curves are the initial density 
profile, (b), (f) Space density of stars initially associated with the primary galaxy; thick curves are the initial density profile, (c), (g) Space 
density of all stars. Lower thick curves are the initial density profile of the primary galaxy, and upper thick curves are the superposition of 
the initial density profiles of the primary and secondary galaxies. Lines of logarithmic slope —1 and —2 arc also shown, (d), (h) Logarithmic 
slope of the surface density profiles of the merger remnants. Thick curves correspond to the initial primary galaxy. 



